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Abstract. We study, using molecular dynamics simulations, the kinetics of shape 
recovery in a model solid undergoing transformations from square to a general rhombic 
lattice, the triangular lattice being included as a special case. We determine the 
necessary and sufficient conditions for such shape recovery in terms of the nature and 
dynamics of transient and localized non-affine zones which inevitably accompany the 
transformation. 

Martensites [I] are classified as reversible (e.g., Nitinol) or irreversible (e.g., steel) 
according to their ability to recover their external shape upon thermal (or stress) cycling. 
What are the conditions under which a martensitic transformation is reversible? A 
recent argument [2] suggests that a necessary condition for reversibility of martensitic 
transformations follows simply from symmetry relations between the parent and product 
phases, viz., reversible martensites are such that the parent and product phases are 
related by a group-subgroup relation, since then a unique parent lattice can be identified 
for every transformed product. On the other hand, irreversibility implies that such an 
identification is impossible or ambiguous. It is clear however that the necessary and 
sufficient conditions for reversibility depend on the dynamics of transformation during 
thermal cycling. 

The conventional approach to the study of the dynamics of martensites assumes 
that the driving force for nucleation can be derived from a non-linear, elastic free-energy 
functional written in terms of a dynamical elastic strain and its derivatives [3]. These 
'strain-only' theories[l] are augmented by the condition of local elastic compatibility, 
restricting the elastic displacements to smooth, single-valued functions. This local 
constraint automatically disallows all configurations with defects and regions of 
plasticity and assumes that the transformation is locally affine at all length and time 
scales. It is known, however, that there is a significant production of dislocations 
during the transformation in irreversible martensites [2]. In order to describe irreversible 
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martensites or the transition from reversible to irreversible behaviour, one needs to go 
beyond strain-only theories of Ref . [I] . 

In this paper, we use a molecular dynamics (MD) simulation of a two dimensional 
(2d) model solid [5j El [7J to (i) check the validity of the symmetry-based criterion (related 
to necessary conditions) and (ii) explore the dynamical conditions for reversibility 
(related to sufficient conditions). This is done by simply changing a single potential 
parameter which allows us to explore both reversible and irreversible martensites within 
the same model. 

Our 2d model solid [SJ El [7J is composed of N particles interacting via a repulsive, 
anisotropic 2-body potential (V2), parametrized by an anisotropy coefficient a, and 
a short-range 3-body potential (V3), of strength t> 3 , which favours local square 
configurations. Specifically, V^r^a) = v 2 (o"o/Vjj) 12 {1 + a cos 2 26^}, where is 
the relative displacement between particles i and j, 6ij the angle between and an 
arbitrary external axis, and V^r^-, r^; V3) = V3 [fijfjk sin 2 + permutations] , with 
Qijk — cos _1 {rjj • *jk/( r ij r jk)}- Energy and length scales are set using v 2 = 1 and o~o = 1 
and the unit of time is a A/m/t> 2 , where m is the particle mass, which for typical values, 
translates to a MD time unit of lps. Decreasing V 3 induces a square (Sqr) to rhombic 
(Rmb) transition while a controls the apex angle of the Rmb phase from tt/6 at a = 0, 
i.e. a triangular (Trg) lattice to 7r/2 for a > 1.2. 

The equilibrium and dynamical properties of structural transitions for this model in 
the constant number - N, area A and temperature T ensemble with fixed external shape 
have been studied in some detail in [H El El- We briefly mention the main conclusions 
below in order to set the stage: 

(i) The transformation from the Sqr to Rmb (or Trg) phase progresses by heterogeneous 
nucleation and growth. 

(ii) The transformation, or order parameter (OP), strain ct (shear for Sqr — > Rmb) 
also produces a non-order parameter (NOP) volumetric strain ey which is slaved 
to Ct- The NOP strain introduces non-local interactions between the transformed 
regions leading to the formation of a twinned microstructure typical of martensites 
for low transformation temperatures. 

(iii) Localized non-affine zones or NAZs are generated at the transformation front where 
the displacement of the atoms cannot be described using affine deformations (viz. 
scaling and shear). In this model solid the NAZs are in the NOP sector. 

(iv) The NAZs in the NOP sector are created when the volumetric stress ay exceed a 
threshold. As soon as the NAZs form, ay tends to decrease thereby screening the 
non-local interactions. At low temperatures, the dynamics of the NAZs in the frame 
of the moving front is slow so that this screening is never total for martensites. At 
high temperatures, on the other hand, the NAZs quickly cover the entire growing 
nucleus destroying the twinned pattern leading to an untwinned, disordered and 
irreversible ferrite. 

(v) Particle trajectories within the NAZs tend to be ordered when martensites are 
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Figure 1. Schematic dynamical phase diagram in T — V3 — a, showing equilibrium 
phases (square (Sqr) and rhombus (Rmb)) by solid lines. The triangular (Trg) solid 
exists for a = V3 = - point P. The black dot marks the location of a tricritical point 
at which the jump in the order parameter vanishes. The dynamical phase martensite 
is formed for quenches at T below the curved surface bounded by the dashed lines. 
Dynamical martensitic phases upon cycling: reversible (RM) and irreversible (IM) 
martensite (red region) are shown. The quench and cycling protocols Qi and Q2 are 
denoted by arrows, (b) Zero temprature energy per particle E/N for p = N/V = 1.05, 
a = and W3 = .408 as a function of the OP strains ei and near the Sqr — > Rmb 
transition showing a metastable Sqr minimum at (0,0) and two degenerate, stable Rmb 
minima, (c) E/N(e' 2 ,e' 3 ) for p = N/V = 1.05, a = and V3 — 3.034 for the reverse 
transformation from the Trg to Sqr phase. The strains are now calculated from the 
Trg phase at U3 = a = 0. 



obtained while they are disordered for the high temperature ferrite. 

In this study we are interested exclusively in shape transformations accompanying 
the formation of martensite so the temperature is set to a low value (T = 0.1) 
throughout. Further we need an isolated solid with stress free boundaries to allow 
for deformations of the shape of the solid as it transforms. We associate every particle 
with a "glue" density, £(r) = 1 for r <r g dropping smoothly to zero at r = r g (chosen 
to be the next-nearest neighbour distance). The embedding energy of particle % in this 
glue is V g = —Kg J2j=i n ^( r «i)' where K g is the cohesive energy The glue causes the 
particles to stick to each other producing a 2d solid whose boundary is self consistently 
determined by the many-body interaction among particles alone [10]. The value of K g is 
tuned so as to maintain the density to be roughly constant across the transformation - 
the area A, on the other hand, may vary at fixed pressure and T. Our MD simulation[9] 
uses a leap-frog Verlet algorithm with a time step of .001 which conserves energy to 1 in 
10 -6 , and a Nose-Hoover thermostat to obtain trajectories of particles in the constant N, 
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Figure 2. (color-online) Configurations of particles as v$ is first reduced and then 
increased taking the system from the Sqr to the Rmb phase and back. The colors 
correspond to the quantity a;, = uiq X^fc sin 2 (A9ijk) where the particles j and k are 
near neighbors of i and wo is chosen so that u> varies from 1 (red) in the Rmb/Trg 
to (blue) in the Sqr phase, (a)-(c) Qi: Configurations at T = 0.1, p = 1.05, 
a = 1.25,JV = 14400. (a) v 3 = 10 (b) v 3 = 0.5 (c) v 3 = 10. (d)-(f) Q 2 : same as 
above with a — with v 3 — 10 (d), w 3 = 0.5 (e) and v 3 = 6 (f). (g)-(i) same as (d)-(f) 
except N = 100, the small size of the system makes the transformation reversible. 



area and T ensemble. This gives us an equilibrium phase diagram inT — v 3 — a[6j [7J E], 
which we display schematically in Fig. 1(a). We start with an equilibrium square solid 
at t — and cycle the control parameters as shown by the arrows ( Qi and Q 2 in Fig. 
1(a)). 

The variation of the T = energy per particle E(e%, e^j/N as a function of the OP 
strains ct = {e2i e z) where e2 is the deviatoric and e% the shear strain is is shown in 
Fig. []](b) and (c). For a general first order p4m — > p2 transition, there are four minima 
apart from the minimum at (0, 0) corresponding to the Sqr phase. For the special case 
of Sqr — > Rmb transition shown in Fig. H](b), C2 = and the four minima collapse 
into two. Note that the product Rmb phases are connected in OP space to an unique 
parent Sqr phase so that the reverse transformation is also unique. If a = 0, (FigH^c)) 
however, decreasing v 3 finally leads us to the Trg structure which has a higher symmetry 
(p6m) than Sqr, there being three symmetry axes. Measuring the strains from the Trg 
lattice, we find that the strain energy now has four minima. The central one at (0,0) 
corresponds to the Trg lattice which is surrounded by the three degenerate Sqr minima 
only one of them being the original parent Sqr phase. 

In our model solid we can change the group- subgroup relation of the parent-product 
by changing a. Keeping a = 1.25 we decrease v% to t; 3 = 0.5, in steps of .5 holding the 
system for 10 4 MD steps at each V3 (Qi in FigJTfa)). The Sqr — > Rmb phase transition 
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Figure 3. (a) Overlap parameter A vs. v% at T = 0.1,N = 14400, and p = 1.05 
for various values of a plotted during the forward (<— , decreasing i> 3 ) and reverse (— ►) 
paths. The meaning of the symbols used are explained in the key to the left, (inset) 
Shape recovery as measured by A cyc / - the difference in values of A after one cycle as 
a function of the size of the solid N for a = .5 K g = 1 (full line) and K g = 2 (dashed 
line). Note that large cohesive energy (K g ) helps shape recovery, (b) NAZs for the 
Rmb — > Sqr reverse transformation for a = 1.25. Darker regions have larger D^. Note 
that the NAZs are confined mainly near the surface. 



at fg w 1. breaks the (Ising) symmetry between the two degenerate Rmb minima. Since 
the presence of the surface breaks translational symmetry nucleation predominately 
proceeds from the surface or the corners. The crystal structure as well as the overall 
shape of the crystallite transform from Sqr to Rmb (Figj2] (a)-(c)). In the reverse path 
when t> 3 is increased again, the shape change reverses in accord with Ref.j2] to Sqr. 
The situation, on the other hand, is quite different when a = shown by the line Q2 
in FigfjTa). The Sqr — > Rmb phase transition occurs at v£ = 1.4 and as V3 is further 
reduced to zero, finally a Trg lattice results. The overall shape of the product crystal 
is not rhombic (Fig^d)-(f)), because the Ising symmetry is not completely broken and 
grain boundaries exists between different degenerate variants of the triangular phase 
forming at different portions of the crystallite, - the grain boundary energy between the 
variants of the Trg product being low. During the reverse transformation, these grains 
tend to transform to different Sqr lattices not necessarily the one which produced them 
in the first place, setting up large internal stresses, which are accomodated initially by 
shape deformations but cause the solid to rupture[12] when excessive. 

The overall shape change of the sample during a quench cycle, may be quantified as 
follows. We first identify boundary atoms by counting the number of nearest neighbors. 
The coordinates of the selected atoms are then used to define the shape function rs(9) 
with the angle 9 measured from the x— axis. Denoting r^(6) as the shape of the initial 
square configuration, we have the overlap A given by, 



where L is the system size and we have taken care to check for configurations related 
to the initial configuration by global translations and rotations. The result is plotted in 
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FigHHa) for two values of a as v% is first reduced from a large initial value (= 10) to 0.5 
and then increased again for a system of N = 14400 particles. As expected, we observe 
that for a = 1.25, the transformation is perfectly reversible. For a = .5 when the 
product resembles the Trg phase, A tends to increase during the reverse transformation 
showing that the shape change becomes irreversible. The value of A at the end of one 
cycle A cyc i measures overall shape recovery at the end of the cycle. In the inset of 
Figl^a) we plot A cyc i for a system of particles at a = .5 as a function of N. In view of 
the conclusions drawn from Ref. [2], it is striking that the shape transformation, while 
irreversible for large N, becomes reversible as N decreases. Increasing the cohesive 
energy K g from 1 to 2 also reduces A cyc [ and aids shape recovery. 

We now demonstrate that all of these observations may be understood from 
the dynamics of NAZs which inevitably accompany the solid - solid structural 
transformation [7J. As in Ref. [Zj we study the dynamics of the NAZs by computing 
the local non-affine parameter using the following procedure. For every particle in the 
Trg configuration, we define a neighbourhood Q, using a cutoff equal to the range of the 
potential (~ 2.5a ). This is compared with that of the same particle in the transformed 
lattice by defining the parameter |13j. 



which needs to be minimized with respect to choices of affine strains e mn . Note that 
the indices mandn = 1,2 and r"(0) and r"(t) are the n th component of the position 
vector of the i th particle in the reference (Sqr) and transformed lattice, respectively. 
The residual value of Dq(t, t) is a measure of non-affineness. 

The transformation connecting the Sqr and Rmb phases is accompanied by NAZs 
as shown in Figj3](b). Firstly, as in [5 J and [7j these NAZs are in the NOP sector, 
being associated with the local volume strain ey. They are mainly localized near the 
surface and rapidly disappear as the crystal transforms from Sqr to the Rmb phase, 
being advected out at the completion of the transformation. Secondly, as in Refs.[5] 
and [7], particles close to the NAZs within the transformed region move ballistically 
and in a coordinated manner. It is these two properties of the NAZs discussed here 
that ultimately renders the Sqr — > Rmb martensitic transformation reversible, in spite 
of significant transient and localized plastic deformation. 

Even the slightest amount of plasticity in the OP sector, on the other hand, would 
make the transformation irreversible. Within our model system, a deep quench to the 
a = 0, v 3 = region produces a Trg solid which is not related to the parent square 
lattice by a group-subgroup relation[3]. During the reverse transformation, therefore, 
there is no unique parent lattice that the system can revert to. This produces non- 
affineness in the OP sector due to a multiplicity of affine paths and destroys reversibility. 
The result from our MD simulations is shown in FigJU^a). It is interesting to note 
that the largest values of D 2 exists along boundaries of isolated patches within which 
the transformation is to a single square lattice. Particles at the boundary of these 
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Figure 4. (a)The non affine parameter Dq (see text) for a Trg lattice during the 
reverse transformation at u 3 = 5, a = 0, p = N/L 2 = 1.05 and T = 0.1. Dark regions 
correspond to large D a within non-affine zones (NAZ)s which surround isolated and 
disjoint regions where Dq is small, (b) Particle trajectories in a NAZ - red square in 
(a). Note that individual trajectories are disordered, (c) The same system as in (a) at 
a later time and ^3 = 7.5. Note that the system fractures along lines with high Dq, 
the crack nucleates at a spot on the surface, shown by the red circle in (a) , where Dq 
is particularly large. 



patches are structurally frustrated and tend to follow separate possible affine paths 
along trajectories which are disordered and "diffusive" despite the temperature being low 
(FigH^b)). Unlike in reversible martensites, the creation of these NAZs, do not reduce 
stress. Indeed, the stress continues to increase being set by the OP strain and eventually 
the crystal fractures (FigH^b)). Examination of the fracture surface shows that the solid 
breaks apart precisely along regions of large D Q so that the NAZs provide seeds for the 
heterogenous nucleation of cracks. In essence, therefore, microstructural reversibility in 
martensites, is related to the nature of the accompanying plastic deformation. 

Is is now easy to see why shape recovery is restored for smaller and stiffer solids. 
Since NAZs are produced when the local stress increases beyond a threshold, the average 
distance between NAZs is determined by the statistics of the stress threshold and of the 
local stress. In general one expects that NAZs are separated by some average lengthscale 
l c determined by the yield strength of the material, with stiffer materials having larger 
l c . It may thus be possible to avoid OP strain induced NAZs completely either by 
reducing the system size or increasing l c by making the material stronger. 

In this paper, we have examined the problem of shape recovery of a solid undergoing 
a structural transition from a Sqr to either Rmb (p4m — > p2) or Trg (p4m — > p6m) 
lattice. We obtain necessary and sumcent conditions for shape recovery by examining 
local regions of plasticity viz. NAZs produced during transformation. In agreement 
with Ref. [2 J we show that a group - subgroup transformation is always reversible since 
it is accompanied by reversible NAZs related to the slaved NOP strains. For group 
- nonsubgoup transformations, in general, shape change is not recoverable due to the 
presence of plasticity in the OP sector, unless the size of the system is smaller than the 
typical distance between NAZs. We believe that our work has relevance to applications 
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of real martensitespQ as well as on the general question of phase ordering dynamics of 
solid state transformations. 
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